Non-exponential London penetration depth in Bai_a;K2;Fe2As2 single crystals 
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We have studied the in- and out-of-plane magnetic penetration depths in the hole- doped iron 
based superconductor Bai_a;Ki,Fe2As2 (Tc ~ 30K). The study was performed on single crystals 
grown from different fluxes and we find that the results are nearly the same. The in-plane London 
penetration depth Xab does not show exponential saturation at low temperature, as would be ex- 
pected from a fully gapped superconductor. Instead, Xab{T) shows a power-law behavior, A oc T" 
(n « 2), down to T ~ 0.02 Tc, similar to the electron doped Ba(Fei_a:Coi,)2As2. The penetration 
depth anisotropy 7a — Xc{T)/Xab{T) increases upon cooling, opposite to the trend observed in the 
anisotropy of the upper critical field, 75 = i/^'^(0)/ff|2 (0). These are universal characteristics of 
both the electron and hole doped 122 systems, suggesting unconventional superconductivity. The 
behavior of the in- plane superfluid density Pab{T) is discussed in light of existing theoretical models 
proposed for the iron pnictides superconductors. 

PACS numbers: 74.25.Nf,74.20.Rp,74.20.Mn 



The discovery of superconductivity {T^ ~ 26 K) in 
fluorine doped LaPeAsO ("1111") [if has generated re- 
markable interest in the community. In a short time, the 
critical temperature was increased by pressure or chemi- 
cal substitution above 55 K, which is significantly larger 
than the highest Tc reported in any s-wave superconduc- 
tor, i.e. MgB2, and comparable to those of the cuprates. 
Later, superconductivity with Tc as high as 38 K was dis- 
covered in Bai_a;Ka;Pe2As2 0]. This BaK-122 compound 
is particularly important because, unlike the cuprates or 
the 1111 iron pnictides, it is not an oxide, downplaying 
the role of oxygen in this type of high temperature super- 
conductors. Moreover, large, high quality single crystals 
of both electron and hole doped 122 pnictides were syn- 
thesized 0, H, [1] , which is essential for drawing reliable 
conclusions regarding their physical properties. 

One key feature for understanding the origin of 
the high critical temperature and the pairing mecha- 
nism in pnictide superconductors is the symmetry of 
the order parameter. Phase diagrams of electron Q 
and hole-doped 122 show that superconductivity 
emerges through doping, suppressing the orthorhom- 
bic/antiferromagnetic (AF) ground state in the parent 
compound. Close proximity to a magnetic state could 
imply the importance of magnetic fluctuations for pair- 
ing and may be reflected in the symmetry of the su- 
perconducting gap. At the same time, band structure 
calculations and ARPES experiments [1, [§| show that 
multiple bands cross the Fermi level, opening the possi- 
bility for multiband superconductivity. Several ARPES 

studies 0, El, m 

on Bai_a;K2;Fe2As2 have found at 
least two different superconducting gaps without nodes 
in the a6-plane. Point contact spectroscopy llSll. spe- 
cific heat jl4'| and microwave penetration depth [l5| data 



also suggest fully gapped superconductivity, and possi- 
bly two gaps. However, ^-SR studies show a linear 
temperature dependence of the superfluid density at low 
temperatures and spin-lattice relaxation rate 1/Ti, from 
As NMR [l3| was found to be proportional to T^. To- 
gether with more recent reversible magnetization mea- 
surements [Si, these results leave open the possibility 
for a nodal gap. Extensive work on the electron doped 
122 via Co substitution (FeCo-122) have revealed 

that at low temperatures the penetration depth does not 
show exponential saturation, but instead it features a 
robust power-law behavior AA(T) oc T", with n being 
between 2 and 2.5, depending on the doping level. 

In this paper, we report measurements of the Lon- 
don penetration depth in the hole doped 122 compound 
Bai_a;K;j,Fe2As2. To account for the effect of sample 
preparation on experimental data and draw objective 
conclusions, single crystals of Bai_:rKa;Fe2As2 with Tc « 
30 K, grown by two different groups, using different 
fluxes, were studied. A total of six samples was measured, 
giving highly reproducible results. We show the data for 
two of them. The crystal labeled A, with dimensions of 
740 X 800 X 70/im'^ was grown out of Sn flux. It comes from 
one of the batches characterized in Ref . [5| ■ Based on el- 
emental analysis from Ref. sample A corresponds to 
optimal K doping with a;=0.45, i.e. Bao.55Ko.45Fe2As2. 
Sample B, with dimensions of 320x830x70 /im'^, was 
grown from self-flux of FeAs and it was characterized in 
Ref. 0. According to Ref. [6|, sample B corresponds to 
a::=0.3, i.e. Bao.7Ko.3Fe2As2. 

The penetration depth was measured using a tunnel 
diode resonator (TDR) technique 2l|, by inserting the 
sample into the inductor of a self-resonating tank circuit 
powered by a tunnel diode. The resonant frequency of 
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FIG. 1: (Color online) Xab{T) for sample A, grown from Sn 
flux and for sample B, grown from FeAs self- flux g]. In- 
set: Temperature dependence of resistivity for each sample, 
normalized to the room temperature value. 

the empty resonator was /o « 14 MHz, with a stabil- 
ity better than 5 ppb/hour. The rf magnetic field pro- 
duced by the inductor was H^/ ^ 10 mOe, much less than 
typical values for the lower critical field Hd ~ 50 Oe in 
iron-arsenides. Therefore, when a superconducting sam- 
ple is placed inside the inductor, the inductance of the 
coil changes due to Meissner screening, which leads to a 
change in the resonant frequency. The relative frequency 
change is directly proportional to the penetration depth, 
A/ = — GAA [22]. The calibration constant G was de- 
termined both by following a calibration procedure de- 
scribed in Refs. [l^, and from in-situ extraction of 
the sample from the coil, with both methods yielding 
identical results. First, we mounted the sample with the 
c-axis along the rf magnetic field {Hrf |j c). In this geom- 
etry, only the in-plane screening currents were induced, 
and therefore A/ oc AAot. Then, the sample was aligned 
with Hrf -L c. For this geometry, screening currents flow 
both in the ab plane and along the c direction and AAc 
could be obtained in our relatively thick samples using 
the numerical model described in Ref . [l^ . 

Figure [T] shows Xab{T) (main figure) and resistivity (in- 
set) for samples A and B. We notice that although they 
correspond to different potassium concentrations, both 
samples have very similar values for the onset Tc- The 
widths of the transition at Tc are significantly different, 
with AT < 3K for sample A and AT < 1.5K for sample 
B. Moreover, despite the fact that optimal doping corre- 
sponds to a; = 0.45 K for both methods, the maximum 
Tc reached by using Sn fiux is about 30 K whereas 
by using FeAs flux it is 38 K A possible explanation 
is that some amount of Sn is incorporated in the crystal 
structure of sample A, as explained in Ref. 

Figure shows the low temperature region of Aab(T) 
for samples A and B. Despite the differences in potas- 
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FIG. 2: (Color online) (a) Low temperature region of \ab{T) 
for Sample A and Sample B, respectively. Inset: AXab vs. 
(T/Tc)^ for hole doped (Sample B from present work) and 
electron doped (a;=0.074 Co from Ref. HI]) FeAs-122 (shifted 
vertically for clarity). Dashed lines represent linear fits, (b) 
Aai)(r) for Sample B and our theoretical calculations based on 
the model from Ref. [23|, for the following values of impurity 
scattering parameters (Fq, F^) (top to bottom): (0,0.064), 
(3, 0.064), (3,0.068), (3,0.060). Inset: The lowest temper- 
ature region, showing good agreement between experiment 
and theory. 



slum content and transition width at Tc, the penetration 
depth behavior is similar below 0.3Tc, exhibiting a power 
law Aofc(T) = Ao + 5T", with the average value of the ex- 
ponent being n « 1.9 ± 0.4. It is important to note that 
this nearly quadratic temperature dependence of the pen- 
etration depth is similar to previous observations in the 
electron-doped 122 [ll,|23|. 

The inset to Fig.[l]shows AAq6 plotted against (T/Tcf 
for sample B and for a Ba(Fei_2:Co2:)2As2 single crystal 
(x=0.074) [l^, corresponding to optimal doping. It can 
be clearly observed that at low temperatures, both traces 
are well fit by straight lines with similar slopes. This very 
similar power-law behavior suggests that there are com- 
mon electrodynamic properties within the 122 pnictide 
superconductors. 
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FIG. 3: (Color online) In-plane Aab(T) and out-of-plane Ac(T) 
penetration depth for Sample B. Inset (a): Variation with 
temperature of 7a = Xc/Xab obtained by dividing the values 
from the main panel, and 7^(T) from Ref. [s^]- Inset (b): 7a 
for hole (current work) and electron doped (Re. [3l|]) FeAs- 
122. 



There are several possible explanations for this ob- 
served behavior. The existence of point nodes on 3D 



parts of the Fermi surface was suggested in Refs. 
A quadratic temperature dependence for the penetration 
depth is also expected for a superconducting gap with 
line nodes in the presence of strong (unitary) impurities, 
which will create an additional quasiparticle density of 
states 24|. It has also been proposed that antiferromag- 



netic spin fluctuations could mediate inter-band pairing 
in the iron pnictides. In this case, the order parameter 
is fully gapped but changes sign between different Fermi 
surface sheets. This situation is referred to as the ex- 
tended symmetry 25|. It was found that within this 



approach either the existence of line nodes between the 
two Fermi surfaces (26j and/or impurity scattering (27| 
can explain a power-law temperature dependence of the 
experimentally observed NMR relaxation rate and pen- 
etration depth. In Fig [H^b) we plot the experimental 
data for sample B together with our calculations consid- 
ering the extended symmetry (A occosfca;-|-cosfcj,) in 
the presence of non-magnetic impurities for several val- 
ues of the intra-band (Fq) and inter-band (F^) scatter- 
ing rates |23]. Details of these calculations are given in 
Ref. 1131 , where this model was used for the Co doped 
122 system. With Aah(O) as a fitting parameter, we were 
able to correctly reproduce the experimental data down 
to the lowest measured temperature, as seen in the inset 
of Fig.HJb). 

However, the fit yields Aab(O) « 600 nm while experi- 
mentally it was found to be about 180 nm (ii,^. Fur- 
ther theoretical and experimental investigations may rec- 
oncile this discrepancy and it is possible that the power- 



law dependence of A(T) results from the effect of strong 
inter-band scattering on the order parameter with ex- 
tended s-wave symmetry. However, the possibility for an 
anisotropic or even nodal superconducting gap remains 
open, particularly given the recent report of linear behav- 
ior of A(T) in LaFePO crystals [s^. We only conclude 
that our results for the hole doped Bai_a;Ka;Fe2As2 sug- 
gest unconventional superconductivity in pnictide super- 
conductors. 

Applying the rf magnetic field both parallel to the 
c-axis and along the a6-plane, we were able to di- 
rectly subtract the change in penetration depth associ- 
ated with screening currents flowing in the c-direction, 
i.e. A Ac. Since we are not aware of independent esti- 
mates of Ac at any temperature, we have used the ap- 
proach of Ref. 3lj. High magnetic field experiments 
on single crystals of Bai_xKa;Fe2As2 [1, [l^l found that 
7^ = -ffc2 (0)/-ffc2(0) « 3-5 near T^. Based on the vahdity 
of the Ginzburg-Landau theory in the vicinity of Tc, we 
consider = ~ 3.5 at T = 0.9 Tc- Figure [3] shows the 
resulting Ac(r) on the same graph with Xab(T). Because 
Ac has a weaker temperature dependence at low temper- 
atures, it extrapolates to a significantly higher value at 
r=0, Ac(0) ~ 1800 nm. Therefore, we obtain a penetra- 
tion depth anisotropy 7a ~ 7 near T = 0, which is about 
twice the value of 7a (Tc). In inset (b) of Fig. [31 we show 
7a (T) from the present work together with the result for 



the electron doped Ba(Feo.926Coo. 074)2 As2 from Ref. [31 1. 
In both cases, 7a decreases with temperature, although 
at different rates, which emphasizes again the similarity 
between the hole and electron doped 122 superconduc- 
tors. 

We now compare the temperature dependence of the 
penetration depth anisotropy 7a (T) with that of the up- 
per critical field 7^(7"). In pnictide superconductors, 
i?c2(0) may reach ~ 100 T, making it difficult to estimate 
7^ over the entire temperature range (3^ . Nevertheless, 
it has been found that 7^ decreases from sa 3.5 near Tc to 
« 1.4 at 0.5Tc (13 and we display both -f^{T) and 7a (T) 
in inset (a) of Fig.[3| A possible explanation for the differ- 
ence in the temperature dependence of the anisotropies 
could be given by the two-gap scenario. In a clean su- 
perconductor at T = 0, 7a (0) = {Va)/{Vc)j where Va and 
Vc are the Fermi velocities along the a and c directions. 
The averaging is performed over the whole Fermi surface. 
However, at T = Tc, 7^(Tc) = {A^vl)/(A'^v^), where A 
is the superconducting gap. For two gaps with different 
values, the main contribution to 7a (Tc) is given by the re- 
gions of the Fermi surface with a larger gap. If the larger 
gap exists on a less anisotropic Fermi surface, it will re- 
sult in 7a (Tc) < 7a (0). This situation is the reverse to 
that of MgB2, where 7a increases with temperature (33|. 
However, unlike MgB2, as we show below, a simple two 
gap s-wave model cannot correctly describe the superfluid 
density, pafc(T) = {Xab{0) / XabiT)f in Bai-.K.FesAss. 

Figure [H shows Pab{T) for Sample B calculated for 
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FIG. 4: (Color online) Blue symbols: PabiT) for Xab{0) = 
180 nm. Orange line shows best fit to the two-gap s-wave 
model. Red symbols: pab{T) for Aab(O) = 600 nm. Black 
lines show calculated Pab{T) for the same values of the scat- 
tering parameters, ro/27r Tco and rT/27r Tco, shown in Fig. [J] 
obtained from the fit to extended s-wave model. Inset (a): 
Zoom of pab{T) for Aai,(0) = 180 nm (symbols) with the 
best fit to two-gap model. Inset (b): Zoom of pabiT) for 
Aa6(0) — 600 nm (symbols) with the best fit to s"*" model. 



Fig. [2l Inset (b) zooms on the low-temperature part and 
indicates that the best fit to s"*" model works down to the 
lowest temperature. 

In conclusion, the hole doped compound, 
Bai_a;Ka;Fe2As2, exhibits a nearly quadratic tem- 
perature variation of the London penetration depth 
down to 0.02Tc, at odds with the exponential behavior 
expected for conventional fully gapped superconductor. 
This is similar to the behavior observed for electron 
doped Ba(Fei_2;Co2:)2As2 [3, Efll and suggests uncon- 
ventional superconductivity in the entire 122 pnictide 
family. As a note, we mention that similar measurements 
of rf penetration depth by the Bristol group also found 
Xab{T) oc at low temperatures |35l] . 
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Aa6(0)=180 nm and 600 nm. For the experimental value 
of Aa6(0)=180 nm (blue circles), Pab{T) has positive cur- 
vature at T > Tc/2, whereas both s- and d-wave calcula- 
tions predict only negative curvature in the entire tem- 
perature range. This may be explained assuming two 
independent s-wave gaps Ai(T) and A2(T), with rela- 
tive contributions e and (1 — e), so that p(T)=ep(Ai)-|- 
(1 — e)p(A2). Orange line shows the best fit to this model. 
Even though the quality of the fit looks acceptable at in- 
termediate temperatures, the model clearly fails to repro- 
duce the low temperature data, as illustrated in the inset 
(a) of Fig.m Moreover, the parameters obtained from the 
fit: Ai(0) « OMknT^, A2(0) « 0.25kBTc and e « 0.8 
are physically unreasonable. Both gaps are smaller than 
the BCS weak coupling limit of 1.76 fcsTc and they 
are much smaller than those previously reported from 
ARPES [13, llll |. These two gaps would reproduce mea- 
sured Tc only if inter-band coupling is strong, in which 
case our fitting model of independent gaps is not appli- 
cable [sl]. We have tried the same fit for Aafc(0)=600 
nm and arrived at the same conclusion. Furthermore, 
theoretical calculations using the extended symmetry 
in the presence of impurities do not reproduce the ex- 
perimental data for Aab(0)=180 nm either. On the other 
hand, if we use Pab{T) calculated with Aafc(O) = 600 nm 
(red circles in Fig. [4|) a good agreement with s+ model 
can be achieved in the entire temperature range. Black 
lines in Fig.|3]show Pab{T) calculated for four sets of scat- 
tering parameters, ro/27r Tco and r7r/27r Tco, shown in 
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